Scaling in thermal convection: A unifying theory 



Siegfried Grossmann 1 and Detlef Lohse 1,2 

1 Fachbereich Physik der Philipps-Universitat Marburg, 

Renthof 6, D- 3 5032 Marburg, Germany 

2 University of Twente, Department of Applied Physics, 
P.O. Box 217, 7500 AE Enschede, The Netherlands 

(February 5, 2008) 

A systematic theory for the scaling of the Nusselt number Nu and of 
the Reynolds number Re in strong Rayleigh-Benard convection is suggested 
and shown to be compatible with recent experiments. It assumes a coher- 
ent large scale convection roll ("wind of turbulence") and is based on the 
dynamical equations both in the bulk and in the boundary layers. Several 
regimes are identified in the Rayleigh number Ra versus Prandtl number Pr 
phase space, defined by whether the boundary layer or the bulk dominates the 
global kinetic and thermal dissipation, respectively. The crossover between 
the regimes is calculated. In the regime which has most frequently been 
studied in experiment (Ra < 10 11 ) the leading terms are 
Re ~ Ra 1 / 2 Pr~ 3 / 4 for Pr < 1 and Nu ~ Ra^Pr' 1 ^ 2 , Re ~ Ra 1 / 2 Pr- 5 / 6 
for Pr > 1. In most measurements these laws are modified by additive cor- 
rections from the neighboring regimes so that the impression of a slightly 
larger (effective) Nu vs Ra scaling exponent can arise. The most impor- 
tant of the neighboring regimes towards large Ra are a regime with scaling 
Nu ~ Ra 1 / 2 Pr 1 ! 2 , Re ~ Ra 1 / 2 Pr~ 1 / 2 for medium Pr ("Kraichnan regime"), 
a regime with scaling Nu ~ Ra 1 ^ 5 Pr 1 ^ 5 , Re ~ i?a 2 / 5 Pr~ 3 / 5 for small Pr, a 
regime with Nu ~ Ra 1 ^, Re ~ i?a 4 ' 9 Pr _2//3 for larger Pr, and a regime with 
scaling Nu ~ Ra?/ 7 Pr' 1 / 7 , Re ~ Ra^Pr^l 7 for even larger Pr. In par- 
ticular, a linear combination of the 1/4 and the 1/3 power laws for Nu with 
Ra, Nu = 0.27 -Ra 1 / 4 + 0.038-Ra 1 / 3 (the prefactors follow from experiment), 
mimicks a 2/7 power law exponent in a regime as large as ten decades. For 
very large Ra the laminar shear boundary layer is speculated to break down 
through nonnormal-nonlinear transition to turbulence and another regime 
emerges. - The presented theory is best summarized in the phase diagram 
figure 0. 
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I. INTRODUCTION 



The early experiments on turbulent Rayleigh-Benard (RB) convection in air cells with 
Prandtl number Pr pa 1, summarized by Davis |l|], showed a power law increase of the 
Nusselt number Nu with the Rayleigh number Ra, namely, Nu ~ Ra 1 with 7 = 1/4. 
[] However, in these early experiments only relatively small Rayleigh numbers Ra < 10 8 
were achieved. Later, when RB experiments with larger Rayleigh numbers and in water 
cells with Pr pa 7 were done, the power law exponent 7 turned out to be larger than 1/4. 
Malleus' elegant theory of marginal stability, resulting in 7 = 1/3, seemed to describe those 
experiments [[|. 

In the late 80's, Libchaber et al.'s experiments done at the University of Chicago on high 
Rayleigh number Rayleigh-Benard convection in a helium gas cell with Prandtl number 
Pr pa 1 revealed new and unexpected scaling for the Nusselt number as a function of 
the Rayleigh number, namely Nu ~ Ra 1 with 7 = 0.282 ± 0.006 [^, |J. The Reynolds 
number Re, characterizing the wind near the walls, i.e., the large eddy mean flow, scaled as 
Re ~ Ra a with a = 0.491 ± 0.002 These results were reproduced and extended in many 
experiments and numerical simulations |5 H23|1 ; for review articles, which also summarize the 



results of earlier experimental, theoretical, and numerical work, we refer to refs. [T0,[T7 



24fl . From all these experiments at first sight it seems that at least the scaling exponent 
7 pa 0.282 ± 0.006 pa 2/7 is very robust. 

Various theories were put forward to account for the scaling of the Nusselt number Nu 
and the Reynolds number Re as functions of the Rayleigh number Ra and the Prandtl 
number Pr. These include the Chicago mixing zone model P| and the Shraiman-Siggia 
theory |25|], both reviewed in refs. [K|[T7],[24[]. The main result of the Chicago model is 



Nu ~ Ra 2 / 7 . The Chica go group did not focus on the Pr dependence as only experiments 



with Pr pa 1 were done ||]. Later, Cioni et al. added the Pr dependence in the spirit of 
the Chicago model and obtained 

Nu ~ Ra 2/7 Pr 2/J , (1) 
Re fluct ^Ra^Pr- 4 ' 7 . (2) 

Here, Re fi uct refers to the velocity fluctuations and not to the large scale mean velocity (often 
denoted as the "wind of turbulence") as Re does. This Prandtl number dependence is only 
expected to hold for Pr < 1 [|17],|24|]. On the other hand, the Shraiman-Siggia model, which 
assumes a turbulent boundary layer (BL) and a thermal boundary layer nested therein (i.e., 
Shraiman and Siggia implicitly assume a large enough Prandtl number), states 

Nu ~ Ra 2/7 Pr- l/7 , (3) 
Re ~ Ra 3 ^ 7 Pr~ 5 ^ 7 with logarithmic corrections. (4) 

For an extension of the Shraiman-Siggia theory to position dependent shear rates see ref. 
26|. The Prandtl number dependence of the Nusselt number resulting from the Chicago 



With the symbol "~" we mean "scales as" throughout the text, not "order of magnitude". The 
prefactors are determined in section IV. 
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model 0] and from the Shraiman-Siggia theory |25j are not contradictory, as eq. (|l|) has been 
suggested for small Pr number fluids and eq. (|^) for the large Pr number case. Indeed, in 
the large Pr-limit Zaleski [24| derives the same Pr-dependence eq. @ as in the Shraiman- 
Siggia theory also from the Chicago model. However, both theories are based on rather 
different assumptions. 

In recent years, the Prandtl number dependence of the Nusselt number has been mea- 
sured and comparison with the two theories [fi],p5| became possible. The first experiments 
were done with water and helium RB cells. However, comparing Nu in water and in helium 
convection only allows for a small variation of the Prandtl number. From such experiments 
a small decrease of Nu with increasing Pr at given Ra was reported |27 

76 ±11 for Pr 



monte et al. 27 measured Nu 



0.7 and Nu 



For Ra = 10 9 Bel- 
48 ±6 for Pr = 6.6. Much 



larger Pr variations are possible if RB convection in mercury or liquid sodium is studied. 
Those experiments with a mercury RB cell (Pr = 0.025) by Rossby p8f , by Takeshita et al. 
15| , and by Cioni et al. [ll,|i~7 and with a liquid sodium RB cell (Pr = 0.005) by Horanyi 



et al. |29| reveal that Nu increases with Pr, which is consistent with the Cioni et al. [17] 
extension of the Chicago mixing zone model. 

However, there seems to be indication that also the Chicago mixing zone model cannot 
account for all phenomena observed in recent experiments: One of the most startling ob- 
servations is that there seems to be a small but significant trend of the scaling exponent 
7 as a function of Pr. For Pr w 5 - 7 (water) one has 7 = 0.28 - 0.293 [|K||23,[3(||31| 
for Ra up to Ra « 10 9 and an even larger 7 w 1/3 for larger Ra ~ 10 9 — 10 11 |[32|| : for 
Pr = 0.7-1 (helium gas) it is 7 = 0.282 ± 0.006 gQ; for Pr = 0.025 (mercury) 7 = 0.247 
|H and 7 = 0.26 ±0.02 Q (for Ra < 10 9 ) have been measured; and for Pr = 0.005 (liquid 
sodium) it is 7 = 0.25 P9"| . Those and further experimental results are summarized in table 



|. All exponents between 0.25 and 0.33 have been measured! For thermal convection in a 
water RB cell (Pr ^5 — 7) with self similarly distributed balls on the top and bottom wall, 
the scaling exponent can even be as large as 7 = 0.45, presumably depending on the ball 
size distribution |33 . 

Next, a breakdown of the 7 ~ 2/7 scaling regime at very large Ra has recently been 
observed, possibly towards a scaling regime Nu ~ Ra 1 / 2 , which has been predicted by 
Kraichnan decades ago |57|| . For Pr = 0.025 Cioni et al. |T7| see the breakdown at Ra 



2 ■ 10 9 (and a startling small window with a local scaling exponent smaller than 2/7 for 
Ra ~ 5 • 10 8 - 2 • 10 9 ) while for Pr = 0.7 - 1.0 Chavanne et al. [0 observe it at Ra ~ 10 u . 
The transition around Ra ~ 10 11 in figure 3 of ref. ||10|| ) showing Nu/ Ra 2 / 7 vs Ra from Wu 
and Libchaber's data P, |3"8"|1 , may already be interpreted as the same breakdown. On the 
other hand, Glazier et al. [[3(J do not observe such a transition. Thus, the experimental 
situation itself is not yet clear. 

All these observations and also the more intuitive rather than equation of motion based 
approach of ref. Q call for a re-examination and extension of the existing scaling theories 
for thermal convection. Of course, a mathematically rigorous derivation of Nu(Ra,Pr) 
and Re(Ra, Pr) is hardly possible. The known rigorous bounds overestimate the measured 
Nusselt numbers by more than one order of magnitude and are only able to give the scaling 
exponent of the Kraichnan regime 7=1/2 [|39| -|4T ] . 



Though in shortage of a strict mathematical derivation, the guideline of the presented 
approach will be the dynamical equations for the velocity field u(x, t), the kinematic pressure 
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Reference 



fluid 


Pr 


Ra range 


7 


SF 6 


1-93 


10 9 - 10 14 


0.30 ±0 


H 2 


5.5 


10 7 - 3 • 10 9 


0.293 


H 2 


6.8 


3 • 10 7 - 4 • 10 9 


0.290 


H 2 


6.5 


10 9 - 2 • 10 11 


1/3 


H 2 


ps 7 


2 • 10 8 - 2 • 10 10 


0.28 


H 2 


ps 7 


2 • 10 8 - 2 • 10 10 


0.28 ±0 


H 2 


ps 7 


8 • 10 7 - 7 • 10 9 


0.281 ±0 


He 


0.8 


4 • 10 5 - 2 • 10 9 


0.280 


He 


0.7-1 


< io 11 


0.282 ±0 


He 


0.6 - 1.2 


4 • 10 7 - 10 12 


0.285 


He 


0.6-0.73 


3 • 10 7 - 10 11 


2/7 


air 


PS 1 


< 10 s 


0.25 


Hg 


0.025 


2 • 10 4 - 5 • 10 5 


0.247 


Hg 


0.025 


10 6 - 10 8 


0.27 


Hg 


0.025 


5 • 10 6 - 5 • 10 8 


0.26 ±0 


Hg 


0.025 


4 • 10 8 - 2 • 10 9 


0.20 


Hg 


0.025 


2 • 10 5 - 8 • 10 10 


0.29 ±0 


Na 


0.005 


< io 6 


0.25 



Ashkenazi & Steinberg 34] 
Garon & Goldstein 



Tanaka & Miyata [31] 
Goldstein & Tokuda [pi 
Qiu & Xia 
Lui & Xia |2 
Shen et al. gj 
Threlfall @ 
Castaing et al. [| 
Wu & Libchaber fl 
Chavanne et al. Jl 
Davis I! 



Rossby g!8 
Takeshita et al. [15| 
Cioni et al. [17 
Cioni et al. f7) 
Glazier et al. fl36l 



Horanyi et al. [29| 



TABLE I. Power law exponents 7 of the power law Nu ~ Ra 1 for various experiments. The 
experiments were done with different aspect ratios, however, no strong dependence of the scaling 
exponent 7 on the aspect ratio is expected (in contrast to the prefactors). 



field p(x,t), and the temperature field 9(x,t), 

d t Ui + UjdjUi = -d^ + vd]ui + (3g5 i3 9, (5) 

d t e + Ujd s e = Kd]e, (6) 

assisted by the appropriate boundary conditions at the bottom wall z = 0, the top wall 
z = L, and the side walls of the cell. Here, g is the gravitational acceleration, j3 the isobaric 
thermal expansion coefficient, v the kinematic viscosity, k the thermal diffusivity, and L the 
height of the RB cell; the temperature difference between top and bottom walls is called A. 

The second feature of our approach introduced in section II is that we try to be as 
systematic as possible. We will be able to identify four different main scaling regimes for Re 
and Nu in the Ra — Pr phase space, depending on whether the BL or the bulk dominates 
the global thermal and kinetic energy dissipation, respectively. Three of the four regimes 
consist of two subregimes, depending on whether the thermal BL or the viscous BL is 
thicker. We also calculate the validity range the scaling laws and make predictions on the 
stability of the different regimes. In section III we compare the power law exponents of the 
theory with experimental data. In section IV we try to adopt the prefactors of the theory to 
some experimental information and compare the resulting prefactors to further experiments. 
Section V contains a summary and conclusions. 
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II. BOUNDARY LAYER VS BULK DOMINANCE OF KINETIC AND 

THERMAL DISSIPATION 



A. Definitions 

The parameter space of RB convection is spanned by the Rayleigh and by the Prandtl 
numbers, 

Ra = Pr = 1. (7) 

KV K 

Our main focus is on the resulting Reynolds and Nusselt numbers, 

where (.) A denotes the average over (any) z-plane. Correspondingly, (.) v used below denotes 
the volume average. U is the mean large scale velocity near the boundaries of the cell. It 
is the remainder of the convection rolls which in the turbulent regime manifests itself as 
coherent large scale convection flow, as first discovered by Krishnamurti and Howard W2 



and later found by various groups [|, [10|, |2C|, [27|, |38|, |43H46f . The existence of this "wind of 
turbulence" is one of the central assumptions of our theory. We consider this to be a 
weak assumption, given the overwhelming experimental evidence. The effect of the wind is 
twofold: (i) In the range between the wind and the cell wall a shear flow boundary layer 
will build up. (ii) The wind stirs the fluid in the bulk. In the presented theory we consider 
the velocity fluctuations in the bulk of the cell only as a consequence of the stirring by the 
large scale roll. Therefore, the Reynolds number Re based on the roll velocity rather than 
the one based on the fluctuations Reluct i n the bulk is taken as the more appropriate one 
to theoretically describe the bulk turbulence. 

B. Decomposition of the energy dissipation 

The starting points of the present theory are the kinetic and thermal dissipation rates 

e u (x,t) = u(d i u j (x,t)) 2 , (9) 
e e (x,t) = K (di6(x,t)) 2 . (10) 

Their global averages (e u (x,t)) v = e u and (ee(x,t)) v = eg obey the following rigorous 



relations, which are easily derivable from the equations of motion, see e.g. [10, 25 



e u = ^(Nu-l)RaPr- 2 , (11) 
A 2 

e d = k—Nu. (12) 

Dissipation takes place both in the bulk of the flow and in the boundary layers. Near the walls 
thermal and kinetic boundary layers of thicknesses Xg and X u develop, which are determined 
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by the thermal diffusivity k and the kinematic viscosity v, respectively, and which are in 
general different, depending on Pr. They are defined on the basis of the temperature and of 
the velocity profiles, respectively. Whenever there exists a thermal shortcut in the bulk due 
to the turbulent convective transport, the width of the thermal boundary layer is connected 
with the Nusselt number by 

X e = l -LNu- 1 . (13) 

The thickness of the kinetic boundary layer can be expressed in terms of the Reynolds 
number, 

X u ~ LRe~ 1/2 . (14) 

Here, we have assumed that there is laminar viscous flow of Blasius type (cf. sects. 39 and 41 
of ref. ||47|| ) in the boundary layer; the lateral extension x of the BL has been identified with 
the height L of the cell, reflecting that the wind organizes in the form of a large scale roll. 
The transition to turbulence in the boundary layers will be considered in subsection II-F. 
Though for large enough Ra the total volume of the BLs is rather small, their contribution 
to the global average dissipation may be considerable, as the velocity and the temperature 
gradients in the BLs are much larger than in the bulk. 

In general, we decompose the globally averaged dissipation rates into their BL and bulk 
contributions, 

^u,BL ^ufiulki 
Eg = £0,BL + Eg^bulk, 

where e U)B l = J < z <\ u + Il-x u <z<l dzp(diUj) 2 /L = u {(diUj(x G BL,t)) 2 ) v is the viscous 
dissipation taking place in the viscous BL, Eg^L = Io<z<x fl Jt ~ Sl-\ b <z<l 

dzn{did) 2 /L = 

k {{did{x G BL,t)) 2 ) v is the thermal dissipation taking place in the thermal BL, e u ^ u ik = 
J Xu<z<L _ Xu dzu(diUj) 2 / L = v ((diUj(x G bulk,t)) 2 ) v is the viscous dissipation taking place 
in the bulk, etc. 

This kind of thinking immediately suggests the existence of four regimes: 

(I) Both e u and Eg are dominated by their BL contributions; 

(II) eg is dominated by e 9)B l and e u is dominated by 

Eu ,bulk} 

(III) e u is dominated by e U)B l and eg is dominated by eg^ik', 

(IV) both e u and eg are bulk dominated. 

For (relatively) small Ra the BLs are most extended, therefore regime I is expected. On 
the other hand, for large Ra the BLs are very thin and we will expect regime IV, provided 
that the volume reduction of the BL is more efficient than the dissipation increase in the 
BLs due to the growing shear rate. Next, for small Pr the viscous BL is smaller than the 
thermal one, X u Xg, and we expect regime II. Finally, for large Pr it is A„ ^> Xg and we 
have regime III. 

A priori it is not clear whether all four regimes can exist. However, after input of some 
experimental information, we will see that they are likely to exist. We will also calculate 
the scaling of the border lines between the different regimes in the Ra — Pr phase space. Of 
course, these lines do not indicate sharp transitions but the range of the change of dominance. 
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C. Estimate of bulk and BL contributions 



The next step is to estimate the various contributions e Uj BL, £e,BL, £u,buiki e e,buik of the 
BL and the bulk dissipation from the dynamical equations (0) and expressing them as 
functions of Nu, Re, Ra, and Pr. 
Bulk contributions: 

We start with the kinetic dissipation. As already outlined above, the theory assumes that 
the bulk fluctuations with typical velocity u fi uct originate from the large scale coherent flow 
with velocity U . If there is no such "wind of turbulence" , the following estimates cease to 
be valid; the Reynolds number Re = UL/v even cannot be defined properly Clearly, this 
will happen in the regime of very large Prandtl numbers where the flow is suppressed by the 
strong viscosity. Therefore, there will be a transition line Pr(Ra) in phase space, defined 
by, say, Re = 50, beyond which the theory no longer holds; we will calculate this line in 
section II-E. Another limit of applicability is in the very small Pr range. Here, k is so large 
that the heat is molecularly conducted, thus Nu = 1. Other possible limitations of the basic 
assumptions of the theory will be discussed in section V. 

The assumption of the large scale velocity U stirring the bulk implies the picture of 
a turbulent energy cascade in the bulk which in turn suggests how to estimate the bulk 
dissipation rates, namely by balancing the dissipation with the large scale convective term 
in the energy equations following from (|5|) and (||), 

e u ,buik = v((diUj(x E bulk,t)) 2 J v ~ — = jjRe 3 (17) 

As argued before, we took as the relevant velocity scale the wind velocity U and not the 
velocity fluctuations Ufi uct because it is U which stirs the fluid in the bulk. This is a key 
assumption of the theory, justified by intuition and by the results. The presented theory 
does not make any statement on the -Ra-scaling of the typical bulk fluctuations Ufi uct . 

We explicitly remark that the findings on the Re dependence of the energy dissipation 
rate in Taylor-Couette flow by Lathrop et al. do not contradict eq. fllTD - Lathrop et al. 
|48| found that e u L 4 z/ _3 i?e -3 still depends on Re even for large Re. However, their result 
refers to the global e u , not to e u f,uik- ~ Possible intermittency corrections are not taken into 
consideration in eq. (|17|) as they are at most small . 



We note that strictly speaking there should be a factor (L — 2X U )/L on the rhs of eq. 
(|T7j), as the average e u ^ u ik = (e u (x G bulk,i)) v refers to the whole volume. However, it is 
assumed that the state is already turbulent enough, i.e., Ra large enough, so that X u <C L. 
The validity of this assumption limits the scaling ranges to be derived. 

The estimate of the thermal bulk dissipation ee^ u ik is slightly more complicated, as the 
velocity field u(x,t) matters in the dynamical eq. (|6]) for the temperature. In particular, it 
matters whether the kinetic BL, characterized by a linear velocity profile, is nested in the 
thermal one or if it is the other way round. 

For the former case (X u < Xg, i.e., small Pr, the thermal boundary layer can be estimated 
in complete analogy to eq. ( |TTD as 



ee,bui k = « ((di0(x e bulk, t)f) y ~ — = K—PrRe. (18) 
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Note that in correspondence with U in eq. ([171) we took the large scale temperature dif- 
ference A in eq. (|IBD , not the typical temperature fluctuations A^ nct . Again, this is a key 
assumption, justified by the later results. 

For the latter case (X u > Xg, i.e., large Pr, we must realize that at the merging of the 
(linear) thermal BL into the thermal bulk the velocity is not U itself, but smaller by a factor 
Xg/X u < 1. Therefore, it is reasonable to assume that UXg/X u is the relevant velocity for 
the estimate of eg^uik, i-e., 

Xg UA 2 A 2 



e<9,fru2fc 



X u L L 2 

BL contributions: 



n—PrRe^Nu' 1 . (19) 



For e U: BL w e follow an idea by Chavanne et al. Jl9| and estimate, using eq. (|14D, 

e u ,BL = v{(d iU3 (x E BL,t))% ~ ^ ~ ~ ^ Re5/2 ' (20) 

Here, U/X u characterizes the order of magnitude of diUj and the factor X u /L accounts for 
the BL fraction of the total volume. Again, this reasoning breaks down when there is no 
large scale "wind of turbulence" . Correspondingly, we estimate 

e 9 ,BL = k ({diOix E BL, t)f) v ~ ■ ^ ~ k^Nu. (21) 

Eqs. flT7p to ([2"T]) express the various dissipation contributions (and thus the total dissipations 
e u and eg, (|I"5| ) and (|T6|)) in terms of Ra, Pr, Re, and Nu. If we insert eqs. (|17|) to (pTf) 
into the rigorous relations ( [TT]) and (P^), we obtain two equations, allowing to express iVu 
and i?e in terms of i?a and Pr. If we only take the dominating contributions e BL or e W(! 
in e u and eg, respectively, the formulae for the four regimes I, II, III, and IV are obtained, 
describing pure scaling instead of superpositions. 

With this idea in mind, the scaling of the thermal boundary layer dissipation fl2T|) , though 
correct, does not give new information. It coincides with the rigorous relation (p!2|). The 
physical reason is that the bulk is considered to provide a thermal shortcut. Therefore, we 
make use of the dynamics in the thermal BL in more detail. We approximate (systematically 
in order 1/Re) eq. @ by the dominant terms (cf. or [17, p5]) 



u x d x e + u z d z e = Kd;e (22) 

in the thermal BL. Both terms on the lhs are of the same order of magnitude as can be 
concluded from the incompressibility condition d x u x + d z u z ~ 0. In the lower subregimes 
with X u < Xg the velocity u x must be estimated by U , in the upper subregimes with X u > Xg 
it is as argued above u x ~ UXg/X u . In addition, d x ~ 1/L and k<9 2 ~ k/A^. Therefore, for 
X u < Xg we finally get 

Nu ~ Re l/2 Pr 1/2 (23) 

and for X u > Xg we have 

Nu ~ Re 1/2 Pr l/3 . (24) 

Eqs. (^) and (^|) replace the correct (but useless) relation eg ~ eg^L which does not add 
new information beyond (|12^). 



8 



D. Four regimes 



We will start with the eg^uik dominated regimes (III and IV). 
Regime IV, e u ~ e uMlk and e e ~ e eMlk (large Ra) 

Depending on whether X u is less or larger than Xg we must use eq. (|18D or eq. (|1J), re- 
spectively, for the ee^u/fc-estimate. The former happens for low Pr, the latter for large Pr. 
Therefore, we will give these two subregimes the index "1" for lower and "u" for upper. At 
what line Pr(Ra) in phase space X u = Xg, the crossover from the X u < Xg to the A„ > Xg 
will occur is not clear a priori. We will later calculate this line X u = Xg with additional 
experimental information. 

In regime IV\ we use ( |T?D for e u in QTTD and for eg in flT2| ) to obtain 

Nu ~ Ra l/2 Pr 1 ' 2 , (25) 
Re ~ Ra 1/2 Pr~ 1/2 . (26) 



We recognize the asymptotic Kraichnan regime [0, just as expected for large Ra when 
both thermal and kinetic energy dissipation are bulk dominated. Note that other lines of 
arguments can also lead to eq. fl25|), see e.g. Kraichnan's work itself |37[], Spiegel [0, or our 
reasoning in section II-F. Therefore, eq. (p5|) seems to be quite robust. The physics of this 
regime is that the dimensional heat current NukA/L is independent of both k and v. 

In regime IV U we substitute as before ( |17D for e u into eq. (|TTD but now (|I9|) instead of 
(P~S|) for eg into eq. (|TJ) to obtain 

Nu ~ .Ra 173 , (27) 
i?e ~ Ra A/9 Pr- 2 ' 3 . (28) 

The A/tt scaling is the one also following from the Malkus theory 0. 
Regime HI, e„ ~ e u , B L and 6g ~ ggj^fc (large Pr) 

Again we have to distinguish between the lower subregime with X u < Xg and the upper 
one III U with X u > X e . For UI X we combine (M) with O) and (|l|) with (03) and get 



iVu ~ Ra 2/3 Pr l/3 , (29) 
i?e ~ Ra 2/3 Pr~ 2/3 . (30) 

This regime will turn out to be small and less important. The more important one is III U ' 
Combine (M) with ([11]) and (0) with QT3) to obtain 



A^m ~ Ra 3/7 Pr- 1/7 , (31) 
i?e ~ Ra^Pr-*' 7 . (32) 

This regime may be observable for large enough Pr when A„ ^> Ag. To our knowledge up 
to date this regime has neither been observed nor predicted. 

Later, we will find hints for this regime III U in form of a subleading correction to describe 
the Chavanne et al. data [|19|]. It would be nice to perform further experiments with large Pr 
to be able to more cleanly identify this postulated regime III U . We note that this regime is 
not in contradiction to Chan's |51| upper estimate Nu < constRa 1 / 3 , holding in the infinite 
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Pr limit (for fixed Ra), and also not in contradiction to Constantin and Doering's [52 
rigorous upper bound Nu < constRa l l 3 {l + log Ra) 2 / 3 , holding in the same Pr — > oo limit. 
The reason is that regime III U is for finite Pr; if Pr — > oo, also Ra — > oo, if one wants to 
stay in regime IH U - 

Regime II, e u ~ e Ujbldk and e e ~ t e , BL (small Pr) 



Regime I If. Combining (TT7]) with (|TTJ) gives together with ( |2"3"| ) 



A^m ~ Ra 1/5 Pr 1/5 , (33) 
i?e ~ Ra 2/5 Pr~ 3/5 . (34) 



This regime should show up for small enough Pr when A u <C Xg. Indeed, Cioni et al. [17 
observed experimental hints for such a regime; also eqs. (^) - (|34D have already been derived 
by them in a similar way |T7]]. A power law Nu ~ i?a 1/5 was already suggested by Roberts 



Regime Because of the two competing conditions eg ~ £g,BL (i- e -, -P r small) and X u > Xg 
(i.e., Pr large) such a subregime can at most be small. It will turn out later that it will 
probably not exist altogether. Nevertheless, for completeness we give the scaling laws, 
resulting from taking now (p4|) and, as before, inserting (|17]) into (|i~l~D, namely 



Nu ~ Ra 1/5 , (35) 
Re ~ Ra^Pr- 2 ' 3 . (36) 



regime 


dominance of 


BLs 


Nu 


Re 


/i 


e u ,BL, £e,BL 


X u < Xg 
X u > Xg 


0.27i?a I / 4 Pr 1 / 8 
OMRa^Pr- 1 / 12 


OWPa^Pr" 3 / 4 
OmdR^^Pr" 5 / 6 


Hi 
(IIu) 


£u,bulki e 6,BL 


X u < Xg 
X u > Xg 


OmRa^^Pr 1 / 5 
(~ Ra 1 / 5 ) 


0.47Pa 2 / 5 Pr~ 3 / 5 
(~ Pa 2 / 5 Pr- 2 / 3 ) 


nil 

I IIu 


£u,BL, eg,bulk 


X u < Xg 
X u > Xg 


6.43 • 10~ 6 i?a 2 / 3 Pr 1 / 3 
3.43 • 10~ 3 Ra^ 7 Pr~^ 7 


5.24 • l(T 4 Pa 2 / 3 Pr- 2 / 3 
6.46 • 10~ 3 Pa 4 / 7 Pr- 6 / 7 


IVi 

IV U 


€u,bulki £0, bulk 


X u < Xg 
X u > Xg 


4.43 • lO^Ra^Pr 1 / 2 
0.038RC1 1 / 3 


QmQRa}l 2 Pr- l l 2 
0.16Pa 4 / 9 Pr" 2 / 3 



TABLE II. The power laws for Nu and Re of the presented theory, including the prefactors 
which are adopted to four pieces of experimental information in section IV. The exact values of the 
prefactors depend also on how the Reynolds number is defined, see the first paragraph of section 
IV. Regime II U is put into brackets as it turns out that it does not exist for this choice of prefactors. 

Regime I, e u ~ e u , BL and e e ~ e 9jBL 

Regime This is the regime of (comparatively) small Ra whose scaling we obtain from 
using (|23"D and substituting ( p0|) for e u in (jlTf ), namely 

Nu ~ Ra^Pr 1 / 8 , (37) 
Re ~ Ra 1/2 Pr- 3 '\ (38) 
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FIG. 1. Phase diagram in the Ra — Pr plane. The power laws and the corresponding prefactors 
(to be determined in section IV) in the respective regimes are summarized in table ||. The tiny 
regime right of regime is regime IIIi. On the dashed line it is X u = Xe- In the shaded regime 
for large Pr it is Re < 50, in the shaded regime for low Pr we have Nu = 1. The dotted line 
indicates the nonnormal-nonlinear onset of turbulence in the BL shear flow discussed in subsection 
II-F. The scaling in regime II', is therefore as in the bulk dominated regime IV\ . - The power laws 



for the boundaries between the different regimes are given in table III 



We argue that this is the regime whose scaling behavior has been observed in almost all 
thermal turbulence experiments |3|-|T|, p|-pT|,p!5|-p0|. |22|. [T6|] , but that in nearly all cases the 
pure scaling behavior ( |37D and (^) has been polluted by sub-dominant contributions from 
the neighboring regimes, as we will elaborate in detail in the next section. 

Remarkably, it is this power law Nu ~ Ra 1 ^ which was the first one suggested [0J and 
which has been well known in the engineering literature for a long time j54|- It also holds 
for two-dimensional convection in the low Prandtl number limit [|55|,|56 . 
Regime I u : The scaling in I u is obtained from equation (|4|) and combining (|20|) with fllTl), 
namely 



Nu ~ Ra^Pr" 1 / 12 , 
Re ~ Ra 1/2 Pr- 5/6 . 



(39) 
(40) 



Note that the Ra dependence is the same as in Ii, but now Nu decreases with increasing Pr. 
This behavior is physically to be expected because due to increasing v the convective heat 
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transport is more and more reduced. And indeed, such a crossover from increase to decrease 
of Nu with Pr has been observed in experiment. It will later give us the opportunity to 
determine the transition line X u = Xg. 




FIG. 2. Nu vs Pr for fixed Ra = 10 6 . The fit is based on all points. The expected power law 
Nu ~ Pr 1 / 8 i n regime l\ is also shown (dashed). For larger Pr > 2 one enters regime I u where 
Nu ~ p^ 1 / 12 (dashed) is expected. 

The scaling of this crossover line can already be determined here from equating Xg ~ 
L/Nu and X u ~ L/^/Re in the respective regimes. We obtain 

Pr^- /u ~ Ra°, (41) 
Pr [h-iiu „ Ra ^ (42) 

Pr nh-in, „ i2a -i/ 2j (43) 

Pr m-iv u „ Ra ~i/^ (44) 

Because the line X u = Xg obeys Pr = const in regime I it either is above or below the 
common corner point of all four regimes. Therefore it can go either through regime III or 
through regime II, but not through both. Thus either regime II U will exist or regime IIIi, 
never both of them. 

We now calculate the scaling of the boundaries between the other different domains in the 
Ra — Pr phase space. The boundary between I and II is obtained by equating a u ,B l ~ e«,6?iZfe, 
those between I and III by equating eg^L ~ ee,buifc, etc. The results are 
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o r h-lh T> n I (AF\\ 

'trans ±LU 'transi \ ) 

'rtrans 1 ' ~ %rL- ( 46 ) 

iiii-iVi p i / 47 \ 

trans riU 'trans , > \ ') 

IIl-IVi /'4Q\ 

'trans ll,U/ transi \^°) 

Ju-IHu ^ D„3 /' Z .Q\ 

trans ±>M transt l^^/ 

'Hu-IVu ^ D„2/3 /r n x 

rims ll,Ul trans- \ uu / 



Note that all these lines indicate the range of smooth crossover in the dominance of either 
the BL or the bulk dissipation. 

The phase diagram in Ra — Pr phase space with the various regimes and crossovers is 
shown in figure [I], anticipating the prefactors of the power laws, whose exponents we have 
evaluated up to now. We will determine the prefactors in eqs. (p3|) - ( |50| ) in section IV from 
four pieces of experimental information. These experimental informations all come from 
experiments with an aspect ratio of the RB cell of the order of 1. We expect the prefactors 
to depend on the aspect ratio; therefore, all prefactors given in this paper only refer to 
aspect ratio order of 1 experiments. 



E. Range of validity of power laws 

What is the range of validity of the power laws summarized in table |T|? For too small 
Reynolds numbers towards larger Pr, say, Re cr u = 50, the distinction between the bulk and 
the boundary layer is no longer meaningful; the bulk will no longer be driven to turbulence 
by a large scale velocity U. Correspondingly, if the Nusselt number approaches 1 because 
of too small Pr, the splitting of eg in eq^bl an d eg^uik becomes meaningless. Finally, for 
Nu = 1, we no longer have thermal convection but pure thermal diffusion. 

Therefore, we impose the restrictions Re < 50 towards large Pr and Nu > 1 towards 
small Pr. The lines Re = 50 and Nu = 1 are included in above phase diagram figure |1|. 
Their analytical forms directly follow from the power laws of table |I[ they are given in table 

m 

Beyond these lines, in the shaded areas in figure [I], the flow is viscosity dominated or 
thermal diffusivity dominated and the proposed power laws for Re and Nu (table ||) no 
longer apply. 



F. Turbulence transition of the laminar boundary layer 

For very large Ra the theory outlined here requires an extension. It is based so far on 
the existence of a laminar boundary layer flow of Blasius type; its thickness therefore obeys 
X u ~ LRe^ 1 / 2 , cf. section 39 of ref. ||47|| . The shear in this boundary layer is determined by 
the large scale velocity U of the thermal roles in the RB cell and the boundary layer width 
X u . We define the corresponding shear Reynolds number as 

Re s hear = ~ ^/lie. (51) 
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FIG. 3. Nu vs Ra for fixed Pr = 1 (characterizing helium) according to the presented 
theory, i.e., equation (|56| ) (solid line). Also shown is the practically indistinguishable result 
Nu = 0.22i?a ' 289 of a linear regression of (^) (dashed line), which mimicks a power law with 
an exponent close to 2/7. 

The key issue now is that the laminar shear BL will become turbulent for large enough 
Re shear- The details of the mechanism of this turbulence transition is still under study but 
it seems to have nonnormal-nonlinear [|57H 63|| features. What however is agreed upon is 



that the shear Reynolds number at which the turbulence sets in depends on the kind and 
strength of the flow distortion. A typical value for the onset is 



Re s hear,turb — 420. (52) 

It will turn out that such high shear Reynolds numbers can only be achieved in regimes II 
and IV; in the regimes I and III the large Reynolds numbers necessary for the breakdown of 
the laminar shear BL are not achieved. With the information from table |TJ we obtain the 
corresponding line in the Ra — Pr parameter space indicating the laminar-turbulence onset 
range, namely 

Pr turb ~ Ra t £ b (53) 
in regimes Hi and IV U (with different prefactors) and 

Prturb ~ Ra] urb (54) 
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in regime IV\. The corresponding prefactors will be calculated later, but for clarity we 
already included the characteristic lines which mark the onset to turbulence in the BL in 
the phase diagram figure [l] as dashed lines. Above, the Rayleigh-Benard rolls are still laminar 
in the boundary layer, below the boundary layer is turbulent. 

What power laws for Re and Nu are to be expected in the regime beyond the turbulence 
transition of the laminar BL? One might argue that the destruction of the BL laminarity 
means that both the kinetic and the thermal dissipation rates scale as in the turbulent bulk. 
This implies that the scaling of Re and Nu should be the same as in the bulk dominated 
regime IV. In the phase diagram we called those regimes II[, IV/, and IV^. 

The same result is obtained by yet another argument: In a turbulent thermal boundary 



layer it holds 19, 47 



^l^ Nulog (^l\. (55) 



K 



With logarithmic precision the typical velocity scale of the fluctuations in the BL is equal 
to the wind velocity U and therefore eq. fl5"5]) implies Nu ~ RePr. On the other hand, it still 



holds e u ~ e u ,buik or NuRaPr 2 ~ Re 3 . From these two relations one immediately obtains 
the power laws (^) and fl26|), i.e., scaling as in regime IV\ for all three primed regimes II[, 
IV/, and IVl 

Still we feel that further study is necessary to obtain reliable insight about the dissipation 
rate scaling in turbulent boundary layers. This might influence the scaling exponents in the 
primed regimes II[ and IV{ U . 



III. COMPARISON WITH EXPERIMENT: SCALING EXPONENTS 

A. Nusselt number 

The first type of results of the theory which we would like to compare with experiments 
are the scaling exponents. First, we focus on the Nusselt number. 

For fixed Ra = 10 6 the Nu number seems to increase up to Pr pa 7, see figure The 
fit to the experimental data between Pr = 0.005 and Pr pa 7 gives Nu ~ p r o.i3±o.o2 j n 
good agreement with the predicted exponent 1/8 in the regime //. We note, however, that 
the suggestion eq. (|l|) by Cioni et al. |I7| is also consistent with experiment in the small Pr 



regime. 

The increase with Pr seems to cease for Pr between 1 and 10. As stated above, other 
experimental data p7| even suggest a decrease of Nu with increasing Pr in that regime. 



This is compatible with and explained by the present theory which gives a (Ra number 
independent) transition from the T regime with Nu ~ Pr 1 / 8 to I u with Nu ~ p r -^/^^ both 
for fixed Ra. 

As there is even a controversy in literature whether Nu for water with Pr = 6.6 or Nu for 
helium gas with Pr = 0.7 is larger, it is hard to say where exactly the transition from T to I u 
takes place. According to Cioni et al. |T7| it is fair to say that within experimental accuracy 



Nu(Pr = 6.6) = Nu(Pr = 0.7). We adopt this point of view and use it to calculate the 
transitional Pr number Pr^~ Iu to be about 2. This experimental information thus defines 
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boundary between 


Prtrans 


h ~ Hi 


Pr trans = 4.3 • lQ- & Ra t H ns 


T TT T 


trans — l.U ■ 1U Ka trans 


Hi - W 


Prtrans = 9.7 ■ 10 10 Ra^ 


Ilk ~ IVi 


Prtrans = 9.1 ■ 10" 12 Ra\ rans 


In ~ 1 1 In 


Prtrans = 5.7 • 10 i?a trans 


1 1 In ~ IV U 


Prtrans = 4.8 • W~ 8 Ra t ^ ns 


Il-Iu 


Pr x = 2.0Ra° trans 


nil - niu 

IV - IV U 


Pr x = 5.3 • 1CP Ra££ 
Pr x = 7.3 • 10 3 22a^£ 


J r (Re=50) 


-P?Van S = 6.7 ■ 10~ 5 Ra^ns 


J u -(Re=50) 


Prtrans = 3.0 ■ lO" 3 ^^ 


7 r (Nu=l) 


Prtrans = 3.5 ■ 10 Ra trans 


/J,-(Nu=l) 


Prtrans = ^-^Rdtrans 


m M -(Re=50) 


Prtrans = 2.9 ■ lO" 5 /^,,, 


II r (BL-turbul.) 


Prtrans = 9.3 ■ 10 ^ R&trans 


IVHBL-turbul.) 


Pr^ns = 3A-W- 1G Ral ans 


IK-(BL-turbuL) 


Prtrans = 2.3 ■ 10 i?<2 trans 



TABLE III. Boundaries between the various regimes I through IV, towards the limiting regimes 
where Nu = 1 (small Pr) and Re = 50 (large Pr), and the nonnormal-nonlinear onset of shear 
turbulence (last three lines, see section II-F). 
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FIG. 4. The Nu number, as it follows from eq. (56, compensated by two different power laws 
Ra 1 ^ (solid, as suggest by the present theory for the low Ra regime) and Ra 2 /' 1 (dashed). The 
second is hardly distinguishable from a straight line, i.e., pure 2/7 scaling. The inset shows the 
local slope following from eq. (|56"|). 



the line A u = \o in the phase space and separates the lower subregime // with X u < Xq from 
the upper subregime I u with with \ u > \q. 

Next, we compare the predicted scaling exponent 7 of Nu vs Ra, which in regime I is 
according to the presented theory the same for the lower and the upper subregime. For small 
Pr (mercury, sodium) and (relatively) small Ra the theoretically obtained value 7 = 1/4 
for the scaling exponent of Nu vs Ra has been measured in several experiments: 7 = 0.247 
281, 7 = 0.26 ±0.02 in [ll, and 7 = 0.25 in E3. 
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Also regime II\ (intermediate Ra, low Pr) with 7 = 1/5 seems to have been observed 
by Cioni et al. Q. 

For larger Pr (helium, water) the measured scaling exponent is larger, 7 « 2/7, see table 
|. Here we argue that this results from the superposition of the scaling in regime I, with 
those in the regimes IV U and To substantiate this, we plot the expected Nu vs Ra 

dependence for Pr = 1, 

0.27Ra 1/4 + 0.038Ra 1/3 , 



Nu 



(56) 



in figure ^. Here we have already made use of the prefactors from table O, which will be 
calculated in the next section. Now we fit eq. (^) with one power law in as large a regime 
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FIG. 5. Re vs Ra for mercury, Pr = 0.025. The solid line shows the theoretical superposition 
0.59i?a 1 / 2 +4.30 J Ra 2 / 5 which is very well fitted by a straight line Re = 3.5i?a 446 (dotted, practically 
indistinguishable from the solid line; the fit interval is Ra = 10 6 to Ra = 4 • 10 9 as in experiment). 
The dashed line presents Cioni et al.'s fit through their data Re oc 72a 0,424 . The prefactors cannot 
be compared because of the different definitions of the Reynolds number in refs. and [fTfl| , see 
section IV, first paragraph. 



as 10 5 < Ra < 10 14 . This fit which is nearly indistinguishable from the superposition 
reads 



Nu = 0.22 J Ra ai8y . (57) 

The power law exponent is very close to 2/7 = 0.286 and definitely consistent with the 
experimental data of table |. Changing the fit regime of course changes the exponent of the 
power law fl5"T|). E.g., for a linear regression in the regime 10 6 < Ra < 10 11 which is typical 
for many experiments one obtains Nu = 0.24i?a 0285 , i.e., an exponent which is even closer 
to 2/7. 

By plotting compensated plots or local slopes d log 10 Nu/d\og w Ra as done in figure £| 



for eq. (|56"D one may be able to get hints that there is no pure power law. Note, that on first 
sight a compensation with Ra 2 / 7 may erroneously even be considered as "better" . 

Chavanne et al. [|l9j find hints for a transition to a regime with a visibly larger scaling 
exponent. According to our theory this could be regime III U or IVi or IV/ or, most likely, 
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a mixture of all of them. No clean scaling exponent could hitherto be determined experi- 
mentally. One reason is that in that regime both Ra and Pr change. We will discuss the 



possible nature of this transition below. Also figure 3 of ref. |10| suggests such a transition. 

We now turn to the large Pr regime. There are very few data for large Pr ^> 1. Recently, 
Ashkenazi and Steinberg |54j] performed convection experiments with SF 6 close to its critical 
point. In these experiments both Ra and Pr change considerably at the same time. To what 
degree the RB convection is still Boussinesq close to the critical point is extensively discussed 
[3§. 

0.22 J Ra°- 3±a03 Pr- a2±004 in 10 9 



in ref. 



Ashkenazi and Steinberg obtain Nu = o.22 J Ra u - 3±au3 Pr- a2±U U4 in 10 9 < Ra < 10 14 and 
1 < Pr < 93. Based on the phase diagram [l] we judge that for these ranges of the Ra and Pr 
numbers we should be in regimes I u and III U . From table |I| we see that the Ra exponent 
of Nu is 1/4 and 3/7, respectively. The measured exponent of 0.3 ± 0.03 in between is 
consistent with this. The Pr exponent of Nu is expected to be in between —1/12 and —1/7, 
slightly smaller (modulus-wise) than the value of 
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-0.2 ± 0.04 reported in 
The results of this subsection clearly demonstrate the importance of superimposing the 
power laws of adjacent phase space regimes. This can really mimick different scaling behav- 
ior, as demonstrated in figures [5] and ^. This characteristic feature holds because the power 
law exponents of the neighboring regimes are rather similar. They will commonly show up 
if data in a crossover range are examined. We emphasize that these crossover ranges ap- 
pear to be rather extended, reaching well into the corresponding regimes, due to the small 
differences of the scaling exponents. 

Therefore, rather than writing pure power laws, one should allow for superpositions. 
Table |TT| suggests 



Nu 



Ra 1/A Pr' 



c IUu Ra 5 / 28 Pr- 15 / 5e + ■ ■ 
c IVu Ra^Pr- 1 / 8 + ■ ■ 



and 



Re ~ Ra^Pr-V* ■ 



1 + 



V 



■j 
j 
J 

"II, 



# a l/4p r 3/8 + . . . 
jRa -l/20p r 3/40 + . . 



(58) 



jRa l/14p r -3/28 + . . . v 
Pa -l/18p r l/12 + . . . 

p r V4 _| 

jRa -l/10p r 3/20 + . . . J 



(59) 



respectively. Here, we have separated the exponents of regime Which of these corrections 
and how many are to be taken depends on the Pr number and on the aspect ratio. Locally, 
i.e., for a limited Ra range, the suggested Nu vs Ra power law exponents 7 = 2/7 |4].|25| 
(cf. figure 2 of ref. [|19|] or figure 3 of |1C|), or 7 = 5/19 can still be considered as 
an appropriate representation of the experimental data. Globally, for larger Ra intervals, 
however, we claim that eq. (1581) is a better description. 



In previous publications Nu, compensated by the expected scaling Ra 2 / 7 , was plotted 
against Ra in a log-log plot, see figure 2 of ref. From that plot one realizes that the 2/7- 
scaling is slightly too steep between Ra = 10 6 and Ra = 10 8 and not steep enough beyond 
the crossover at Ra = 10 11 . The analogously compensated plot with the expected scaling 
(|3~7|) is shown in figure || Now (for Pr w 1) one obtains a horizontal line up to Ra ~ 10 9 , 
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showing that eq. fl37|) nicely agrees with the data. However, beyond Ra ~ 10 9 one observes 
deviations. We suggest that these corrections originate from the different scaling in the 
neighboring regime III U . The reason that it is regime III U (and not regime IV U ) is that 
in the Chavanne et al. experiments the large Ra measurements also have large Pr; the 
trajectory in control parameter space Ra — Pr is not a straight line. At Ra = 10 10 one 
typically has Pr « f , but at Ra = 10 14 Chavannve et al. typically have Pr »s f — 20. The 
power law exponent 5/28, following from table 0, is consistent with the experimental data 
for large Ra, see figure 0. 



For the mercury data of Cioni et al. [17| we do not have such a complication as Pr = 0.025 



is roughly constant for all chosen Ra. As plotted in figure ^ lower curve, we observe a 
straight line up to about 2 • 10 s and then a decay, signaling contributions from regime Hi. 
The power law exponent —1/20 of the correction term in eq. fl5"8|) is consistent with the data 
shown in figure [| 

A more stringent way to test the superpositions of type (|58"D is to make a linear plot 
Nu/iRa^Pr 1 / 8 ) vs Ra 5 / 28 Pr~ 15 / 5e or vs Ra l l l2 Pr~ 1 ' 8 or vs Ra^Pr 3 ^, etc, depending on 
which neighboring regime the corrections originate. This is done in figure 0, assuming, as 
argued above, that the most relevant corrections originate from regime III U . If the theory 
is correct, the data points must fall on a straight line. Indeed, they do so with satisfying 
precision. 

Note that this kind of linear plot is very sensitive to what combinations of Ra and Pr 
are chosen as x and y axes. E.g., plotting Nu/iRa^^Pr 1 / 8 ) vs Ra^Pr 3 / 8 (the subleading 
correction characterizing regime JVj) does not lead to a straight line at all, see the inset of 
figure [7[ Clearly the variable Ra 5 ^ 28 Pr~ 15 ^ 56 on the abszissa is superior, adding confirmation 
that the Chavanne et al. [ITJJ large Ra experiments represent the physics of regime II I u . 



B. Reynolds number 



We now consider the experimental values for the scaling exponents a of the Reynolds 
number vs the Rayleigh number. For Pr m 7 (water) Xin et al. [(SJ find a = 0.50 ± 0.01 
and Qiu and Xia |2(J find a = 0.50 ± 0.02. Both experiments where done in the Ra interval 
between 2 • 10 8 and 2 • 10 10 , i.e., in regime I u , where exactly this power law exponent a = 1/2 
is expected. For Pr « 1 both Castaing et al. Q and Chavanne et al. [19] find a = 0.49 
for all Ra which suggests that possibly the regimes l\ and IVi are seen where this value 
is predicted, or also regime where the exponent is only slightly higher (4/7). For 

Pr = 0.025 (mercury) Cioni et al. [17] find a = 0.424 from a fit to all available Ra. This 
value is in between the derived values a = 1/2 in regime l\ and a = 2/5 in regime Hi. 

We compare this experimental finding Re oc Ra 0A24 (based on a fit to the data in the 
range up to Ra = 4 ■ 10 9 P7| ) with the superposition according to table [IT] 



Re = 0.59Ra 1/2 + 4.30i?a 2/5 , 



(60) 



cf. fig. |5|. In this relatively short Ra interval the theoretical superposition 
hardly distinguishable from its straight line fit 



is again 



Re = 3.5Ra 



0.446 



(61) 
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whose exponent reasonably well agrees with the measured one. Q 

Moreover, also the theoretically obtained Pr number dependence of Re very nicely 
agrees with available experimental information: Chavanne et al. [19| did experiments 
with (slightly) varying Pr. They then plotted RePr 072 vs Ra and obtained the law 
RePr - 72 = 0.0374i?a 1 / 2 . The exponent 0.72 of the Prandtl number was determined by 
minimizing the scattering of points around a straight line in the log-log plot. It very well 
agrees with the calculated Pr scaling exponent 3/4 in eq. ( p8[) (regime //). A possible rea- 
son for the slight difference between theory and data fit is that part of the experimentally 
realized Ra and Pr already belong to regimes III U and IV U , where according to eqs. (|26|) 
and fl28f ) the expected Pr scaling exponent is 6/7 and 2/3, respectively. However, the devi- 
ation is clearly within the experimental uncertainty. In figure |8] we replot the experimental 
Reynolds number data of Chavanne et al. fli~9" |. In fig. |a we show Re/ Pr~ 3 ^ vs Ra. The fit 



gives a i?a-exponent 0.492 ± 0.002 in very good agreement with the theory's exponent 1/2. 
In fig. we display Re/Ra 1 ^ 2 vs Pr. The data fit results in a Pr-exponent — 0.77± 0.01, 
also in excellent agreement with the theoretical expectation which is —3/4 in 1\ and —5/6 
in I u . 

The large Ra values in both figures will turn out to belong already to the regimes III U 
and IV U , where the Pr-scaling is similar. The available range is too small to perform a more 
detailed comparison. 

The only large Pr data available are again those by Ashkenazi and Steinberg |54]]. They 
obtain Re = 2.6Ra 0A3±0 - 02 Pr-°- 75±0 - 02 in 10 12 < Ra < 3 ■ 10 14 and 27 < Pr < 190. The 
expected Ra exponent of Re is between 1/2 and 4/7, distinctly larger than the measured 
one of 0.43 ±0.02. Similarly, the theoretically expected Pr exponent of Re is between —5/6 
and —6/7, also larger than the measured exponent —0.75 ± 0.02. We have no explanation. 



IV. PREFACTORS 



A. Experimental input to determine the prefactors 

To obtain the values of the prefactors in the power laws within the presented scaling 
theory, we need further input from experiment. But what data to choose? As pointed out in 
the introduction, a huge variety of data is around, often disagreeing with each other even in 
the scaling exponents, not to speak of prefactors. Those often vary by as much as 50% from 
experiment to experiment, even for cells with the same aspect ratio. Another reason which 
makes the adoption to one experiment and the comparison to others difficult is that different 
definitions are used for the Reynolds numbers. E.g., Cioni et al. [[HJ define the Reynolds 



number as Re = 4L 2 f p /i>, where f p is a distinguished frequency at the small frequency 
edge in the temperature spectrum. Chavanne et al. |T3 define Re = UodL/v, where uj is a 



2 The absolute values of the Reynolds numbers cannot be compared here, as the Reynolds number 
definitions in refs. |]l^| and [|l9| cannot be transferred into each other, see the first paragraph in 
section IV. 
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FIG. 6. Compensated Nu vs Ra data for Pr from about 1 through about 20 (helium, upper, 
data taken from ref. Jn|; the higher Ra experiments also have higher Pr-number) and Pr = 0.025 
(mercury, lower, taken from ref. |l7fl ). Also shown are the calculated exponents in the large Ra 
regimes. We also drew the theoretical curve Nu = QMRa^^Pr- 1 / 12 + 3.43 • 10~ 3 Ra 3 / 7 Pr' 1 / 1 
with fixed Pr = 3 to demonstrate that it roughly describes the data. If the expected Pr number 
dependence is considered, the agreement becomes even better. 



typical frequency in the cross correlation spectrum of two temperature signals, measured at 
a vertical distance of d = 2.3 mm and 2 cm off the axis of the cell. 

In spite of these difficulties, we decided to calculate the prefactors of the suggested scaling 
laws by employing the following choice for the input information as a reasonable, realistic 
example. Our reasons are first, to be able to draw a phase diagram with more or less realistic 
values. Second, to stress the importance of the prefactors. But we caution the reader that 
our input choice is somehow arbitrary; other possibilities can equally well be rationalized, 
sometimes shifting the various regime boundaries considerably. 

Above, as an input from experiment, we had already chosen Pr\ = 2 as the Prandtl 
number for which Nu is maximal (for fixed Ra = 10 6 , cf. fig. 0). m addition, we will use 
the following experimental information: 

1. The observed transition Rayleigh number for the transition from regime Ii to regime 
IIIi, Ra trans = 10 11 at Pr = 1 0; 
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FIG. 7. Same data as in the previous figure, but now in a linear plot Nu/(Ra 1 ^Pr 1 ' 8 ) vs 
Pa 5//28 Pr~ 15 / 56 , revealing the quality of the superposition eq. (|5S|), high Pr, regime III U . The 
linear fit (straight line) gives Nu/iRa^^Pr 1 / 8 ) = 0.24 + 3.3- 10- 3 J Ra 5 / 28 Pr" 15 / 56 . The data points 
are taken from ref. with kind permission by the authors. Only data points with Ra > 10 6 
are considered. Note that in this plot Pr varies as much as 0.6 < Pr < 100. - The inset shows 
iVn/(i?a 1 / 4 Pr 1 / 8 ) vs Pa^Pr 3 / 8 ; this variable had to be used if regime IVi contributes the most 
relevant correction. The data do not fall on a straight line. A similar failure results with the regime 
/Vu-compensated variable Ra l l 12 Pr~ 1 ^ . 



2. the observed transition Rayleigh number for the transition from regime J/ to regime 



III, Ratrans = 4.5 • 10 8 at Pr = 0.025 117 

3. the experimental values, taken from ref. []19[, for the Reynolds and the Nusselt number 
at the middle point (Pom, P^m) in the phase diagram figure |l|; 

4. and the prefactor 0.0372 of the scaling law Re = 0.0372.Ra 1 / 2 Pr~ 3 / 4 measured in 
regime Ii pll. 



Information (1) specifies the prefactor of the rhs of eq. (which we call c/ ; _/// ; ) to be 
c^-j/jj = Pr 'irans I[ I ' R a trans = 1-0 ■ 10 22 . Information (2) gives the prefactor of the rhs of eq. 
(|45|) (which we call c^-i^) to be cj-,_// ; = Plans' I 'Rcklans — 4.3 • 10~ 8 . The two curves (f45D 
and (Eft) cross at 
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-0.5 0.0 0.5 1.0 1.5 
lg Pr 



FIG. 8. The figure presents data for the Reynolds number as a function of Ra and Pr from ref. 
|l9|| , with very kind permission of the authors. 

a) i?e/Pr~ 3 / 4 vs Ra. The expected slope is 1/2, the linear regression fit (solid line) gives 
0.492 ± 0.002. 

b) Re/Ra 1 / 2 vs Pr. The expected slope is —3/4, the linear regression fit (dashed line) gives 
—0.77 ± 0.01. The agreement of the prefactor is also excellent. According to theory, following 
table it is log 10 (i?e/i?a 1//2 ) = —1.432 — (3/4) log 10 Pr (solid line, hardly distinguishable from the 
dashed one); the fit value for the prefactor from linear regression is —1.413 ± 0.005. 
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(Ra M , Pr M ) = (1.03 • 10 11 , 0.94). (62) 

This middle point (RaM, PrM) is denned by the conditions € u< bl = £u,buik — e «/2 and 
£e,BL = £e,buik = e e/2- Equation (|62| ) specifies the prefactors of the rhs of eqs. ( |4"?| ) and ( ^8|) 
to be cnit-ivt = 9.1 ■ 10~ 12 and c Hl - IVl = 9.7 ■ 10 10 , respectively. We see that pP^ Iu = 2 > 
Ptm — 0.94 so that there is a regime IIIi, but no regime II U , as already anticipated. The 
line Pr^~ Iu = 2 hits the boundary between / and III at 

(Ra M ,,Pr M ,) = (7.1- 10 10 , 2.0), (63) 

which fixes the prefactors of eqs. (H^) and ([43]) to be 5.7 ■ 10~ 33 and 5.3 • 10 5 , respectively. 



Correspondingly, one obtains the prefactors of eqs. floTf ) and (f44|) to be 4.8 ■ 10~ 8 and 7.3 ■ 10 3 , 
respectively. These are the data on which the phase diagram figure [l] is based; they are 
summarized in table [TIT]. 

Apart from regime II Ii all regimes turned out to have at least one decade of extension 
both in Ra and in Pr and should therefore in principle be visible. However, we should always 
expect one or more subleading corrections. Regime III\ will clearly not be detectable. 

We again stress how dependent this phase diagram drawn in figure [I] is on the choice of 
experimental information. E.g., if we had adopted Glazier et al.'s |36| point of view that 



there is no transition towards a steeper Ra dependence of Nu at least up to Ra = 8 ■ 10 10 , 
regimes II and IV would have shifted further to the right or would even not exist at all. 
But as shown in figure |3|, apparent smooth scaling behavior does not allow to exclude a 
transition. 

Making use now of the experimental information (3) and (4) we can calculate the pref- 
actors in the power laws for Nu and Re. From figs. 2 and 3 of Chavanne et al.'s work [19| 



we can extract the Reynolds and Nusselt numbers at the middle point (RaM, Ptm) of the 
phase diagram which touches all four regimes, namely 

Re M = 1.20 ■ 10 4 , Nu M = 2.78 ■ 10 2 , (64) 

which is information (3) above. The definition of the middle point (RaM, P^m), i-e., the 
conditions £ u ,bl = tu,buik = e«/2 and ee,BL = £e,buik = ee/2, allows to calculate the prefactors 
c Eu ,buik, c ee ,buik, and c e ^ B L on the rhs of eqs. (|17|), (|18[), and (g0|), respectively. One obtains 



Nu M Ra M _ 

z - rr M ne M 



Num 

Ce 9 ,buik = 00 p = 0.0123, (66) 



2Pr M Re M 
Nu m Rclm 



1028. (67) 



Finally, the prefactor cn u on the rhs of relation ( |2"3"D is adopted to Chavanne et al.'s [0| 
experimentally determined prefactor (see figure 3 of that paper) in the relation (^), i.e., 
Re = 0.0372 J Ra 1 / 2 Pr~ 3 / 4 valid throughout regime // (information (4) above). In that regime 



e u ,BL = e u ; with eqs. (JTT]) , (|20|), and (|67|) we get 
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c Nu = (3.72 • 10" 2 ) 2 • c eu>BL = 1.42. 



(68) 



The prefactors referring to the upper halve of the phase diagram are calculated from the 
matching conditions for Re and Nu on the \ u = \q line. With eqs. (|65|) - flBlf ) and the 
matching conditions now all prefactors of the power laws in the four different regimes are 
determined. We have summarized all these power laws in table |I[ 

Note that from the condition X u = \q for Pr = 2 the prefactor of eq. ( |H] ) is also 
obtained. It is 0.25, i.e., \ u = 0.25-LRe™ 1 / 2 . Also the prefactors to the Nu = 1- and the 
Re < 50-borders of validity and the crossover lines to turbulence of the laminar BL flow 
automatically follow, and are included into table |I[ 



B. Comparison of the evaluated prefactors to experiment 

We now would like to compare the absolute agreement of the power laws summarized 
in table [II], whose prefactors result from an adoption to above four pieces of experimental 
information, with further experimental data. We first focus on the Chavanne et al.'s IfTDJ RB 



measurements in helium gas. From the previous section we know that in these experiments 
due to the large Pr at large Ra it is mainly regime III U which causes additional contributions 
to regime l\. Therefore, we have plotted the superposition Nu = 0.33i?a 1 / 4 Pr -1 / 12 + 3.43 • 
10" 3 Ra^Pr' 1 ' 1 with Pr = 3 into fi gure 0. We again stress that in the experiments Pr is 
not constant at all. Nevertheless, the data are satisfactorily described. Note that the solid 
curve in figure |6| is no fit! 

A much better way to check whether the obtained prefactors of the theory agree with 
the measured ones is to do a linear regression of the straight line as offered in figure [7]. 
Such a straight line fit gives Nu/iRa^Pr 1 ^) = 0.24 + 3.3 ■ 10~ 3 Pa 5 / 28 Pr- 15 / 56 . The found 
prefactors are in good agreement with the expectation Nu/ (Pa^Pr 1 / 8 ) = 0.27 + 3.43 • 
10~ 3 Ra 5 / 28 Pr~ 15 / 56 from table |TJ. 

We also compare the theoretical prefactors of the Reynolds number scaling in regime 
Ii with experiment. As the theoretical prefactors have been adopted to the experimental 
R e /p r -3/4 vs Ra power law, it only makes sense to check the prefactors in Re/Ra 1 ! 2 vs Pr. 
From table |lj the expected slope is —3/4 and the expected prefactor 0.037. Linear regression 
gives a slope of -0.77 ± 0.01 and a prefactor of io- L413±0 005 = 0.0386, cf. figure §. 

Next, we compare the experimental Nu vs Ra scaling for mercury with Pr = 0.025 with 



theory. The measured relations Nu = (0.140±0.005)Pa a26±a02 pi, Nu = 0.147Pa a257 p 



and Nu = 0.155Pa a27 ]15[ are all in reasonable agreement with the regime l\ expectation 
Nu = 0.17Pa 1 / 4 from table 0. The same holds for a comparison in regime I If. The reported 
experimental fit is Nu = 0.44 ± O.OlSPa 20 ^ 02 |17|], theory gives Nu = 0.46i?a 1/5 , again, 



remarkable agreement of both the power law exponent and the prefactor. Remember that 
the only experimental input from this experiment into the theory is Rains' = 4-5 • 10 8 . 

Let us also check the prefactors of Nu as a function of Pr for fixed Ra = 10 6 . A 
power law fit to all available experimental data points |TT|,|2"9"| included in figure gives 



Nu = (7.8 ± o.5)Pr°' 13±0 ' 02 which is in agreement with the theoretical expectation Nu = 
0.27i?a 1/4 Pr 1 / 8 = S.SPr 1 / 8 in regime Leaving out the data points for water (Pr = 
7) which strictly speaking already belongs to regime I u gives a slightly larger power law 
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exponent and a slightly larger prefactor, Nu = 8.7Pr 016 . Both exponent and prefactor are 
consistent with the theoretical expectation. 

C. Widths of the boundary layers 

As stated above, the present theory also gives the absolute widths of the thermal and 
the laminar viscous boundary layers, 

\ e = 0.5LNu~ l , (69) 

X u = 0.25LRe- 1/2 . (70) 

The results for the widths of the BLs for Pr = 0.025 (mercury) and Pr = 7.0 (water) 
are shown in figure |9|. In both cases there are three regimes involved, namely, Ii, Hi, and 
IVi for Pr = 0.025 and I u , lll u , and IV U for Pr = 7.0. As expected, for the larger Pr 
numbers the thermal boundary layer is always nested in the viscous one which agrees with 
the experimental observations |[27|| . For the lower Pr numbers it is the other way round. 



If the laminar BL becomes turbulent, the Blasius estimate X u ~ LRe~ 1 ^ 2 for its width 
must be replaced by the thickness of the turbulent BL. To give an idea about this length 
scale, we calculate the width y$ of the viscous sublayer of the turbulent BL within the 
Prandtl theory ||47|| , applied to Couette flow. In the large Re limit it holds [f47|. |65 



^ = 1 . 38 !?f^£>. (71) 

Here, k = 0.4 is the experimentally known van Karman constant. We have included yo/L 
in figure |9] for the relevant large Ra. The turbulent viscous sublayer is thiner than the 
laminar BL. One also notes that for the larger Prandtl number Pr = 7 (water) Shraiman 



and Siggia's assumption |25[] of the thermal boundary layer being nested in the turbulent 
one is just fulfilled. For lower Pr this is not the case any more. 

Many experiments justify the identification of the thermal BL width with the inverse 
Nusselt number, \q = 0.5LNu~ l . For a detailed discussion we refer to the review articles or 
to Belmonte et al. [E7I] or to the more recent work by Lui and Xia [22~|. 



The situation is more complicated for the width of the kinetic BL X u . Its measurement is 
experimentally difficult. Moreover, experimental results on X u seem to exist only for regime 
I. 

Belmonte et al. |27j tried to measure X u in an indirect way, namely, through the detection 
of a spectral cutoff frequency in gas convection which in water convection (at one Ra) is 
found to have a similar dependence on the height z in the RB cell as the velocity profile U (z). 
For Pr « 1 they found X u « const in the regime 2- 10 7 < Ra < 2-10 9 and X u ~ LRa~ OA4±om 
in 2 • 10 9 < Ra < 10 11 . We have no idea about the origin of the measured scaling exponents. 



More recently, Xin et al. [JE], Xin and Xia JBBJ, and Qiu and Xia ptj measured the 
thickness of the kinetic BL in a water cell in a more direct way. They define X u as the 
distance from the wall at which the extrapolation of the linear part of the velocity profile 
U(z) equals the maximum velocity U = max z ll(z), the velocity of the large scale wind. In 
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FIG. 9. Widths of the boundary layers X e = 0.5LNu~ l (solid) and X u = 0.25LRe~ 1 / 2 (dashed) 
for mercury (Pr = 0.025, upper) and water (Pr = 7.0, lower). The dotted lines show the thicknesses 



yo of the viscous sublayers of a turbulent boundary layers, calculated according to eq. (|7l|) . Such 
a width is expected beyond the nonnormal-nonlinear transition to turbulence of the laminar shear 
BL (cf. section II-F). 



the interval 2 • 10 8 < Ra < 10 1U they find X u ~ LRa 
and bottom kinetic BL [f46[ |6(: 



,-0.16±0.02 



for the thickness of the top 



and X u ~ LRa~ 026±om for the thickness of the kinetic BLs 
at the side walls |2(J. The first exponent (for the top and the bottom plates) is different 
from the value of this theory X u ~ LRa~ 1 ^. We can only speculate on the origin of this 
discrepancy. Perhaps, if X u is defined as the distance of the velocity maximum to the wall, 
the i?a-scaling would be different. The power law for the thickness of the kinetic BLs at the 
side walls, however, is well in agreement with the expectation X u ~ LRa~ l l A . 

In any case, the experimentally found very weak dependence of A n on Ra supports the 
assumed laminar nature of the kinetic BL. If the width of the BL were identified with the 
width uq of the viscous sublayer of a turbulent BL, one would expect a stronger dependence on 
i?e, namely yo ~ L\og(k 2 Re)/{k 2 Re) |65[|, i.e., when neglecting logarithmic corrections, 
one would have yo ~ L/ Re ~ LRa~ l l 2 . 
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D. Experimental evidence for the turbulence onset in the BL 



According to the presented theory with the chosen prefactors the breakdown of laminarity 
in the shear BL happens at Ra tur b ~ 10 16 for Pr = 1 and Ra tur b ~ 10 14 for Pr = 0.025. 



These values are calculated from table |T|. Hitherto, there are no experiments for these 
regimes. 

However, one may want to argue that the transition to a turbulent shear BL may already 
occur earlier, be it because of a different aspect ratio, or, in view of more recent work |67|,|68|. 
because the critical Reynolds number (which we had assumed to be 420, cf. eq. (0)) i s 
smaller, or because of a different choice of the experimental input information to which the 
prefactors of the theory are adopted. E.g., if one assumes a laminar shear layer with width 
X u = l.72L/\/ r Re as suggested by Landau and Lifshitz (ref. f47j, section 39) for the (related) 
case of a flat plate shear flow, the transition to turbulence in the shear BL already occurs 
at Raturb = 10 13 (for Pr = 4). We note that this is just at about that Ra number where a 
marked transition in the (thermal) dissipative spectral power has been measured by a probe 
placed in the BL J?|]. This transition was towards a weaker increase with Ra. 

In the context of this section we also interpret the above mentioned recent experiment 
by Ciliberto et al. |33| in which the boundary layers are disturbed by constructing a rough 
bottom plate with a mean roughness comparable to the thermal boundary layer thickness. 
The experiments are performed in water. According to the theory of this paper one would 
expect larger bulk contributions to both the thermal and the kinetic dissipation and therefore 
an earlier onset of regime HI' U and IV^. Indeed, experimentally the increase of Nu with Ra 



is found to be much steeper. For the experiment described by Ciliberto et al. [33] the data 
can be fitted to power laws Nu ~ Ra 35 or Nu ~ Ra 0A5 , depending on the features of the 
rough bottom and upper plate. 



V. SUMMARY AND CONCLUSIONS 

We summarize the central ingredients of the theory presented in this paper: The scaling 
laws for the Nusselt number and the Reynolds number are based on the decomposition of 
the global thermal and kinetic energy dissipation rates into their BL and bulk contributions. 
These in turn are estimated from the dynamical equations, taking the wind U as the relevant 
velocity in the heat conduction cell. The resulting estimates are inserted into the rigorous 
relations eqs. (pTTj) and fll2|) for the global kinetic and thermal energy dissipation, respectively. 
Four regimes arise, depending on whether the bulk or the BL contributions dominate the 
two global dissipations. Each of the four regimes in principle divides into two subregimes, 
depending on whether the thermal BL (of width \g) or the kinetic BL (of width X u ) is more 
extended. 

In addition to these main regimes there is a range for very large Pr numbers in which 
the wind Reynolds number is < 50; here the whole flow is viscosity dominated, and the 
theory looses its applicability. There also is the range of very small Pr numbers in which 
the Nu goes down to Nu = 1, and again the theory no longer holds. Finally, for large Ra 
the laminar kinetic BL becomes turbulent. Beyond turbulence onset we feel the flow is bulk 
dominated. 
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All scaling exponents follow from this theory. If one in addition introduces only four 
pieces of experimental information, also all the prefactors can be determined. Therefore 
the theory has predictive power not only for the power law exponents but also for the 
prefactors. These, however, depend on the chosen experimental information input. To nail 
the prefactors more input information for various aspect ratios is necessary. 

The phase diagram of the theory, the main result of this work, is shown in figure |T]. The 
power laws with the prefactors based on the chosen experimental information are summarized 
in table O, the power laws of the boundaries between the different regimes in table |T|. 

A detailed comparison of the theoretical power law exponents and the prefactors with 
the experimental data gives reasonable and encouraging agreement. We emphasize that to 
accurately account for the dependences of Nu and Re on Ra and Pr single power laws are 
often not sufficient, as additive corrections from neighboring regimes can be considerable. 
This can be viewed as one of the main insights obtained in this paper. A particularly striking 
example is that Nu = 0.27 Ra 1 ^ + 0.038Ra 1/3 mimicks a 2/7 power law scaling over at least 
nine order of magnitude in Ra, see figure |3|. 

The theory also offers a possible explanation why a transition to a steeper increase of 
Nu vs Ra is seen in the Chavanne et al. data |19||, but not in the Chicago group data 



la? 05 0] • ^ ma y be t na ^ i n t ne large Ra Chicago experiments the Pr number was smaller 
than in the Chavanne et al. |T^] measurements. Then for the Chicago data one had to 
expect a transition towards regime IV U where the Ra scaling exponent is 1/3 and thus, 
as demonstrated, in superposition with the leading 1/4 exponent, is indistinguishable from 
a 2/7 scaling. In the Chavanne et al. |19j experiments, on the other hand, one has the 
transition to the large Pr number regime III U where the Ra scaling exponent is 3/7 which 
is much better distinguishable from 1/4. Whether this explanation is true remains to be 
seen. 

Finally, we want to stress and discuss one of the basic assumptions of the theory, namely, 
that a large scale "wind of turbulence" exists, defining the Reynolds number Re = UL/u, 
creating a shear BL and stirring the turbulence in the bulk. Clearly, this assumptions breaks 
down in the shaded area in figure [I] beyond the line Re = 50 where the flow is viscosity 
dominated. But even below this line we do not exclude that the convection rolls break 
down and that the heat is exclusively transported by the fluctuations. E.g., for a water cell 
(Pr w 7) Tanaka and Miyata |3l] do not note the wind of turbulence, in contrast to Zocchi 
et al. E3J , who do observe the wind of turbulence in their experiment with only a slightly 



lower aspect ratio. Also all latest experiments with various aspect ratios doubtlessly detect 
the wind of turbulence [f2U], ^7], [RJ whose existence we therefore consider as a weak - and 
in particular controllable - assumption. 

The present theory does not make any statement how the heat is transported from the 
bottom to the top, i.e., whether it is mainly convective transport or mainly transport through 



plumes ||27|| . Both processes may contribute, as both create thermal and viscous dissipation. 

For even larger Prandtl numbers Pr ^> 7, the spontaneous formation of a wind of tur- 
bulence may seem more and more unlikely. To initiate such a wind so that the presented 
theory can be applied and results can be compared we suggest to slightly tilt the RB cell, 
thus breaking the symmetry and creating a preferred direction for the wind of turbulence. 
If this wind of turbulence can be created, we are confident that the suggested theory holds. 
Presently, we do not have any experimental information which contradicts the theory. How- 
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ever, for further verification more experiments will be valuable. 

Clearly, this paper will not bring the turbulent RB heat conduction problem to an end. 
But we hope that the presented new approach will restimulate the discussion. 
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